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摘 要 : 黄土 高 原 流域 地 抠 系统 的 地 和 狐 演 化 特征 十 分 复杂 , 尚 有 诸多 科学 问题 有 待 进一步 深入 研 


究 。 以 往 研究 大 多 集中 在 流域 地 狐 演 化 的 侵蚀 和 发 育 特 征 等 菜 一 方 


,缺乏 从 流域 地 狐 系 统 及 其 


势能 信息 炳 的 视角 深入 剖析 野外 多 岩 土 层 黄土 小 流域 地 组 演化 特征 的 研究 。 为 此 ,基于 系统 论 的 
观点 和 和 方法 ,构建 多 岩 土 层 黄土 小 流域 地 狐 系 统 及 其 势能 信息 粹 的 数学 模型 ,并 以 辛 店 沟 小 流域 
为 例 , 对 其 地 魏 演 化 特征 进行 研究 。 结 果 表明 :(1) 构建 的 野外 多 岩 土 层 黄土 小 流域 地 魏 系 统 的 概 
念 模型 及 其 势能 信息 六 的 数学 模型 能 够 有 效 对 辛 店 沟 小 流域 进行 数值 模拟 ,(2) 以 黄土 侵蚀 作用 
为 主 的 笠 店 沟 小 流域 从 2000 一 2019 年 的 地 狐 演 化 过 程 是 其 势能 信息 痫 的 炉 减 过 程 和 黄土 地 魏 不 
断 侵蚀 的 过 程 。(3) 辛 店 沟 小 流域 的 势能 信息 炉 能 较 好 地 反映 该 小 流域 的 地 氏 演 化 阶段 和 地 狐 侵 


蚀 过 程 。 
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尽管 早 在 20 世 纪 50 一 60 年 代 就 有 诸多 地 貌 学 
者 充分 认识 到 把 一 般 系统 论 应 用 到 地 貌 学 研究 中 
的 重要 性 ,并 在 把 该 理论 引入 地 貌 学 研究 中 做 出 了 
诸多 有 益 尝 试 " ,但 直到 20 世 纪 80 年 代 初 系统 论 
才 在 地 貌 学 研究 中 得 到 较 广泛 应 用 。 近 年 来 ,从 
系统 论 观 点 出 发 ,对 流域 地 貌 进 行 研究 已 得 到 了 广 
泛 认可 。 励 强 等 "提出 了 关于 流域 地 貌 的 相关 数学 
模型 。 马 新 中 等 中 指出 流域 地 貌 系统 是 一 个 与 外 界 
不 断 进行 物质 和 能 量 交换 的 典型 开放 自然 系统 。 

炉 的 物理 意义 最 早 是 用 于 表征 热力 学 系统 中 


mann 和 Higashi 等 ”对 信息 炉 及 其 测量 方法 开展 
了 进一步 研究 。 艾 南山 中 基于 Davis 和 Strahler 地 貌 
发 育 模型 及 信息 炉 的 概念 ,提出 了 地 貌 信息 炉 的 概 
念 及 其 计算 公式 。 岳 天 祥 等 ”在 地 貌 信 息 炉 的 基 
础 上 进一步 提出 地 貌 信息 超 烂 的 概念 。 由 于 地 貌 
PERES SHO AC fa, PS Hh iE h 
沟 沟 头 活跃 度 和 泥石流 危险 性 评价 等 方面 得 到 了 
BES IZ IAL"). ARSE SW FA TB] JE 
DEM 的 信息 容量 变化 研究 “”。 姜 琼 ” 基 于 系统 
论 和 信息 炉 的 原理 提出 了 势能 信息 炉 的 概念 及 其 


体系 混乱 的 程度 。Leopold 等 中 最早 提 出 了 地 貌 炳 
的 概念 ,但 由 于 地 貌 箭 难以 进行 定量 计算 ,因此 并 
没有 得 到 更 加 广泛 和 深入 的 研究 。 近 年 来 , 随 着 数 
字 高 程 模型 (DEM ) 技 术 的 快速 发 展 ,Zhao 等 "通过 
类 比 热 力学 箭 ,提出 了 流域 地 貌 焙 的 概念 及 其 数学 
模型 ,并 成 功 应 用 于 室内 人 工 降 十 条 件 下 均 质 黄土 
流域 地 貌 演化 过 程 研 究 。 

Shannon 最早 提 出 了 信息 炉 的 概念 。Neu- 
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在 均 质 黄 土 流域 地 貌 系统 中 的 数学 模型 ,并 得 出 了 
势能 信息 炉 在 一 定 程 度 上 可 以 反映 流域 地 貌 系统 
内 部 的 物质 侵蚀 和 堆积 情况 等 初步 结论 。 赵 卫 东 
等 3 基于 势能 信息 炉 对 黄土 小 流域 沟谷 网 络 演 化 
的 变化 特征 进行 研究 ,结果 表明 该 流域 势能 信息 炉 
的 变化 特征 和 沟谷 网 络 的 有 序 化 程度 高 度 一 致 。 
但 上 述 结论 仅 基 于 室内 人 工 降雨 模拟 实验 的 均 质 
黄土 地 貌 演 化 过 程 得 出 ,是否 同 样 适用 于 时 外 黄土 
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高 原 真实 小 流域 有 待 进一步 深入 研究。 为 此 ,本 文 
基于 系统 论 的 观点 和 方法 ,构建 多 岩 土 层 黄土 小 流 
域 地 貌 系 统 及 其 势能 信息 炉 的 数学 模型 ,并 以 黄土 
高 原 辛 店 沟 小 流域 为 典型 样 区 ,对 鞋 店 沟 小 流域 地 
貌 系统 的 势能 信息 依 及 其 焙 变 规律 进行 研究 ,并 据 
此 探索 该 流域 的 地 貌 演 化 特征 ,为 进一步 揭示 注 店 
沟 小 流域 地 貌 演化 机 制 和 土壤 侵蚀 机 理 提供 重要 
TAR STE 


1 研究 区 概况 


本 文 以 黄土 高 原 辛 店 沟 小 流域 (图 1) 为 研究 
区 ,开展 多 和 岩 土 层 黄土 小 流域 地 貌 演化 特征 研究 。 
该 流域 位 于 陕西 省 榆林 市 绥 德 县 无 定 河 左岸 的 辛 店 
沟 , 地 理 范围 为 110°16’~110°17'E、37°29'~37°31'N， 
属于 黄土 丘陵 沟 黎 区 第 一 副 区 ,地 表 侵 蚀 较 严重 ”。 
研究 区 气候 干燥 , 属 半 干 旱 区 温带 大 陆 性 气候 ,降水 
的 年 际 变化 大 ,多 集中 在 汛期 的 6 一 9 月 ”i。 土壤 质 
地 玻 松 ,以 黄 绵 土 为 主 , 主 要 分 布 于 地 形 平缓 或 植 
被 较 好 区 域 **。 


2 数据 与 方法 
2.1 数据 来 源 与 处 理 


为 研究 辛 店 沟 小 流域 地 貌 演化 过 程 , 本 文 收集 
了 该 流域 不 同时 期 的 5 期 DEM 数据 ,如 表 1 所 示 。 


(a) 绥 德 县 


号 广 店 淘 小 流域 
边界 


O 绥 德 县 边 


图 1 辛 店 沟 小 流域 地 理 位 置 


利用 ArcGIS 软件 将 上 述 5 期 DEM 数据 统一 投影 至 
WGS_1984_UTM 投影 坐标 系 , 并 将 所 有 数据 重 采样 
为 90 m 分 辨 率 。 

为 查 明 辛 店 沟 小 流域 内 部 岩 土 层 的 空间 分 布 
及 其 密度 差异 ,作者 于 2021 年 8 月 对 该 小 流域 开展 
了 野外 实地 调查 和 宕 土 密度 的 取样 测试 (图 2)。 该 
采样 点 位 于 辛 店 沟 小 流域 主 沟 道 的 中 下 游 ,从 该 图 
中 可 以 看 出 ,流域 下 部 有 明显 的 岩 士 分 层 情况 , 垂 
直方 向 从 下 往 上 依次 为 砂岩 、 页 岩 和 黄土 。 在 宕 土 
取样 点 1_9( 表 2) 测 得 砂岩 和 页 岩 的 密度 分 别 为 
2.80 gcm 和 2.21 gcm”。 整 个 小 流域 内 部 各 取样 
点 黄土 密度 如 表 2 所 示 , 从 表 中 可 以 看 出 ,整个 研究 
区 内 部 黄土 密度 差异 较 小 。 辛 店 沟 小 流域 黄土 密 
度 取 表 格 中 27 个 取样 点 黄土 密度 的 平均 值 ,为 1.32 
gcm’o 
2.2 多 岩 土 层 黄 土 小 流域 地 貌 系统 概念 模型 构建 

黄土 高 原 小 流域 地 貌 系统 是 一 个 与 外 界 不 断 
进行 物质 和 能 量 交换 的 典型 开放 自然 系统 ”。 一 个 
完整 的 流域 地 貌 系 统一 般 拥 有 四 周 、 上 部 和 底部 边 
界 。 为 构建 多 和 涯 土 层 黄土 小 流域 地 貌 系统 的 概念 
模型 ,首先 必须 确定 该 系统 的 边界 。 本 文 以 辛 店 沟 
小 流域 地 貌 系统 为 例 , 构 建 多 岩 土 层 小 流域 地 貌 系 
统 的 概念 模型 。 根 据 辛 店 沟 小 流域 野外 实地 调查 
结果 ,该 流域 地 貌 系统 的 最 底部 为 砂岩 层 , 上 禾 页 
岩层 ,顶部 为 黄土 屋 。 因 此 ,本 文 把 辛 店 沟 小 流域 


(b) 辛 店 沟 小 流域 À 
图 例 
口 辛 店 沟 小 流域 
高 程 /m 
0 0.25 km 1064.030 
上 一 一 EE 818376 


Fig. 1 Geographic location of Xindiangou small watershed 
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表 1 2000 一 2019 年 辛 店 沟 小 流域 DEM 数据 
Tab.1 DEM data of Xindiangou small watershed from 2000 to 2019 


DEM 数 据 年 份 分 辨 率 /m 数据 来 源 
第 1 期 2000 90.0 SRTM V3(https://www.usgs.gov/) 
第 2 期 2007 12.5 ALOS (https://search.asf.alaska.edu/#/ ) 
第 3 期 2009 30.0 ASTER GDEM (http://www.gscloud.cn/search ) 
第 4 期 2015 90.0 TanDEM-X (http://tandemx—90m.dIr.de/ ) 
第 5 期 2019 8.0 91 卫 图 (https://www.91wemap.com/) 


WÈ : DEM 为 数字 高 程 模型 。 


(a) 野外 涯 土 调查 


(b) 野外 黄土 密度 测量 取样 


图 2 野外 宕 土 调查 和 黄土 密度 测量 取样 


Fig.2 Field geotechnical investigation and sampling of loess density measurement 


地 貌 系统 的 上 部 边界 确定 为 整个 辛 店 沟 小 流域 范 
围 内 的 地 形 表面 ,系统 四 周 的 侧 向 边界 为 沿 小 流域 
分 水 岭 构建 的 垂直 于 水 平面 的 垂直 面 , 系统 底部 边 
界 为 该 小 流域 出 口 附近 高 程 最 低 点 所 在 的 水 平 
面 。 该 水 平面 也 是 辛 店 沟 小 流域 地 貌 系统 的 势能 
计算 基准 面 。 根 据 上 述 岩 土 层 分 布 及 系统 边界 , 构 
建 出 辛 店 沟 小 流域 的 流域 地 貌 系统 的 概念 模型 ,如 
图 3 所 示 。 

2.3 SALE RRMA R A SHAS 
型 构建 

2.3.1 多 岩 土 层 流 域 地 貌 系统 的 剖 分 ”流域 地 貌 系 
统 的 势能 是 指 以 流域 底部 边界 为 势能 计算 基准 面 ， 
计算 出 的 整个 流域 地 貌 系 统 边 界 范围 内 的 岩 土 体 
所 具有 的 总 势能 。 为 了 计算 出 整个 多 岩 土 层 流域 
地 貌 系统 的 势能 和 势能 信息 信 , 需 要 对 整个 流域 地 
貌 系统 进行 平面 剖 分 和 垂 向 分 屋 。 首 先 ,利用 DEM 
正方 格 网 剖 分 技术 ,对 流域 地 貌 系统 进行 平面 剖 
分 ,把 整个 流域 在 平面 上 前 分 为 mxm 个 正方 格 网 单 


元 格 。 训 分 单元 格 的 大 小 和 该 流域 DEM 数据 的 单 
元 格 一 致 。 由 于 势能 与 岩 土 体 的 密度 紧密 相关 , 因 
此 ,为 了 准确 计算 具有 不 同 密度 岩 土 层 的 势能 , 需 
要 根据 涯 土 层 密度 的 不 同 ,对 整个 流域 地 貌 系统 在 
垂 向 上 进行 分 层 。 至 此 ,完成 了 整个 多 兰 士 层 流域 
WARS =A SY AE ,使 其 形成 一 个 包含 mxm 
个 基本 单元 体 的 三 维 立 体 训 分 体 ,每 个 单元 体 在 垂 
向 上 可 以 包含 不 同 密度 的 岩层 或 士 层 。 通 过 计算 
出 每 一 个 单元 体 的 势能 并 进行 累加 就 可 以 计算 出 
整个 流域 地 貌 系统 的 总 势能 。 根 据 上 述 训 分 方法 ， 
考虑 到 辛 店 沟 小 流域 中 出 露 的 基 岩 和 黄土 密度 相 
对 均匀 ,无 需 对 岩 土 层 内 部 进行 进一步 垂 向 训 分 。 
因此 ,该 系统 在 垂 向 上 可 放 分 为 砂岩 页 岩 和 黄土 3 
个 分 层 。 

2.3.2 ELAYE FEAE BI” 
F20194F HET yatta RRA AE E EES 
并 成 功 构建 了 均 质 黄土 小 流域 地 貌 系 统 的 势能 信 
息 炳 数学 模型 。 该 模型 被 用 于 室内 人 工 降雨 条 件 
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R2 辛 店 沟 小 流域 黄土 密度 
Tab. 2 Loess density in Xindiangou small watershed 
序号 点 号 高 程 /m 采样 深度 /cm 黄土 密度 /g: cm 
1 1_1 889 26 1.22 
2 1.2 981 35 1.20 
3 1 3 983 35 1.20 
4 14 918 28 1.30 
5 1_5 929 25 1.28 
6 1_6 857 19 1.28 
T 1-7 970 21 1.42 
8 1_8 1015 24 1.25 
9 1_9 824 20 1.37 
10 1_10 781 10 1.52 
11 1-11 888 13 1.38 
12 1-12 899 16 1.33 
13 2_1 856 12 1.28 
14 2_2 869 15 1.35 
15 2.3 938 17 1.33 
16 2_4 921 13 1.42 
17 2_5 879 20 1.42 
18 2_6 944 24 1.33 
19 2_7 867 19 1.33 
20 2_ 8 931 14 1.23 
21 2_9 955 14 1.37 
22 2_10 948 19 1.33 
23 2_11 926 12 1.30 
24 2_12 919 18 1.35 
25 2_13 955 17 1.32 
26 2_14 968 19 1.35 
27 2-15 969 24 1.30 


下 的 均 质 黄土 小 流域 地 貌 系统 的 地 貌 演 化 研究 , 取 
得 了 较 好 的 研究 成 果 。 但 该 模型 不 适用 于 非 均 质 
多 岩 土 层 流域 地 貌 系统 。 势 能 (6,) 的 计算 公式 
如 下 : 


E,= seed he (1) 


RP: g 为 重力 加 速度 (ms”); p 为 室内 均 质 黄 土 
REg cm”); d HAIER Emm); hh 为 每 个 计算 
单元 体 相 对 于 流域 地 貌 系统 出 口 处 最 低 高 程 的 高 
FE (mm). 

HH , ASCE FAIRS A EE it HAY) 
PP TCA A ABE EB , FERDA F RTE 
Be ASAT IMMA SY AGE CAL 

首先 ,根据 野外 实地 调查 测量 的 基 岩 倾向 、 倾 
角 和 厚度 等 相关 信息 ,计算 出 该 多 岩 土 层 流 域 地 貌 


系统 内 部 第 ; 行 第 7 列 放 分 单元 体 的 基 岩 高 度 ( H, ) 
(起 算 点 为 流域 地 貌 系统 的 底部 边界 ) ,计算 公式 
如 下 : 


H.,= 
sina X sinBY + cosa X sinBX -(A’Y, + ABX, + ACZ,) H 
cos B i 
4=snawxsnB B=cosaXsinB C=cosB 
(2) 


式 中 : a 为 流域 地 貌 系统 中 基 岩 的 倾向 (°); B 为 基 
AAH C); Xo, Yo 为 基 岩 取样 点 的 经 度 、 纬 度 (?°); 
Zo 为 基 岩 取样 点 的 高 程 (m) ;本 为 小 流域 出 口 处 最 
低 点 的 高 程 m);X 了 为 第 ; 行 第 17 列 训 分 单元 格 左 
下 角 的 经 纬度 坐标 (°)。 

其 次 ,多 岩 土 层 流域 地 貌 系统 的 势能 (6,) ,计算 
公式 如 下 : 
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边界 范围 。 

图 3 辛 店 沟 小 流域 地 貌 系统 的 概念 模型 
Fig.3 Conceptual model of the geomorphic system of 
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b= SY sent = Hija- ,) (3) 
式 中 : g 为 重力 加 速度 (ms ”); p, 为 岩层 密度 (g: 

3) ;qd 为 单元 格 尺寸 (m) ; H, 为 第 i 行 第 j 列 剖 分 
单元 格 的 第 k 层 岩层 高 度 (m); 有 ,;, | 为 第 i 行 第 j 列 
剖 分 单元 格 的 第 k-1 层 岩层 高 度 (m) ;为 剖 分 单元 
格 中 宕 土 层 层 数 ;m 为 单元 格 行 数 ;n 为 单元 格 
列 数 。 

姜 琼 3 的 研究 表明 ,流域 地 貌 系统 的 势能 分 级 
BS ARE IC AK , 随 着 势能 分 级 数 
不 断 增加 ， 势能 信息 灶 最 终 会 逐步 趋 于 一 个 稳定 
值 , 此 时 的 势能 分 级 数 为 最 佳 势能 分 级 数 。 利 用 式 
(1) 对 每 个 单元 体 的 势能 值 进行 计算 ,并 根据 势能 
值 大 小 的 不 同 将 其 等 间距 分 为 看 干 势能 等 级 。 最 
JE OT) it Ht FA Be AN a A HE ER CH, ) 进 行 
计算 ” ,计算 公式 如 下 : 

n= piap (4) 


SUP : p() AN TCL i HR ARE ERTE AIERT 
有 的 势能 等 级 中 出 现 的 概率 ;m 为 单元 格 行 数 ;n 为 
单元 格 列 数 。 

利用 上 述 多 宕 土 层 流域 地 貌 系统 的 势能 信息 
烂 数学 模型 ,以 辛 店 沟 小 流域 地 貌 系 统 为 例 , 按 照 
以 下 步骤 计算 出 其 势能 信息 粹 。 首 先 , 根 据 野外 实 
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流域 侧 向 边界 ” ”地 调查 结果 , 辛 店 沟 小 流域 底部 砂岩 层 和 页 岩层 倾 
流域 上 部 边界 > 向 和 倾角 相同 ,岩层 厚度 均匀 ,日 其 厚度 已 经 测 
出 。 利 用 5 期 辛 店 沟 小 流域 DEM 数 据 和 式 (2) 及 式 

a (5) 可 计算 出 每 期 DEM 的 流域 地 貌 系统 势能 (E)。 


m a 1 p? H} + gpa HH, + A + 
>» T 


t 


AFL gp ab, HH Nb, +H, +H) 


SUH sg Ay HD RE m.s”); pi 为 砂岩 密度 (g: 
em”); p, HAARE cm); p 为 黄土 密度 (g- 
cm ) ;d 为 单元 格 尺 寸 (m); ;为 第 i 行 第 7 列 单元 
体 的 砂岩 厚度 (m); 5; 为 第 i 行 第 j 列 单元 体 的 顶部 
界面 相对 于 流域 地 貌 系统 出 口 最 低 点 的 高 度 (m); 
H, 为 页 岩 厚度 (m) ; m 为 单元 格 行 数 ;n 为 单元 格 
列 数 。 

其 次 ,利用 式 (4) 计 算出 辛 店 沟 小 流域 地 貌 系 
统 内 部 的 势能 信息 粹 。 为 了 确定 该 流域 地 貌 系统 
的 最 佳 势能 分 级 数 ,分 别 计 算出 每 期 数据 在 不 同 势 
能 分 级 数 下 对 应 的 势能 信息 箭 , 结 果 如 图 4 所 示 。 
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图 4 辛 店 沟 小 流域 势能 信息 丧 随 分 级 数 变化 
Fig.4 Change of potential energy information entropy 
varying with classification number in Xindiangou small 


watershed 


从 图 4 中 可 以 看 出 , 当 势 能 分 级 数 不 断 增 大 时 ， 
势能 信息 信 与 势能 分 级 数 的 联系 逐渐 前 弱 , 最 后 趋 
于 一 个 常数 。 因 此 ,为 了 准确 计算 出 该 流域 地 貌 系 
统 内 部 的 势能 信息 依 , 需 要 找到 一 个 最 佳 分 级 数 ， 
即 对 势能 信息 nee 吉 果 几乎 没有 影响 的 分 级 
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数 。 本 文 根 据 文献 “了 有 取 势能 信息 炉 的 变化 率 为 负 
万 分 之 一 时 对 应 的 分 级 数 为 最 佳 分 级 数 , 经 计算 最 
佳 分 级 数 为 300 级 。 即 当 分 级 数 为 500 级 时 ,5 期 
DEM 数 据 对 应 的 流域 地 貌 系统 的 势能 信息 粹 的 变 
化 率 很 小 ,可 以 忽略 不 计 , 此 时 可 获得 小 流域 地 貌 
系统 唯一 稳定 的 势能 信息 焙 。 因 此 ,本 文 的 计算 均 
以 500 级 作为 最 佳 势能 分 级 数 。 


3 结果 与 分 析 


3.1 Bae fa SAH eh iE 

SHE a AM E RH TRAR Tit Bh Be PH 
ETE 22 BAr EY Te LEI, IAJ eA 
力 系统 中 分 子 混乱 程度 的 特性 类 似 ”。 此 外 ,势能 
言 息 炉 的 变化 能 反映 流 域 地 貌 系统 内 部 势能 的 变 
化 ,这 与 热力 学 炉 的 炉 变 情况 代表 热力 系统 中 分 子 
运动 能 量变 化 的 特性 类 似 ””。 因 此 ,势能 信息 烂 的 
变化 能 在 很 大 程度 上 反映 流域 地 貌 系统 内 部 黄土 
侵蚀 特征 和 势能 空间 分 布 特征 。 姜 琼 ” 利 用 势能 
言 息 炉 对 室内 人 工 降 雨 条 件 下 均 质 黄 土 小 流域 (以 
下 简称 “室内 小 流域 ”) 的 黄土 侵蚀 特征 进行 了 初步 
研究 ,结果 表明 势能 信息 箭 的 依 减 过 程 能 够 有 效 反 
映 黄土 小 流域 的 侵蚀 过 程 。 为 进一步 验证 该 结论 
是 否 适 用 于 野外 真实 多 兰 土 层 黄土 小 流域 ,以 注 店 
沟 小 流域 地 貌 系统 为 例 , 人 研究 其 势能 信息 信 与 黄土 
侵蚀 之 间 的 关系 。 


(a) SABE fies MMi SY) A PE 
2.1, ”一 - FRE —— Yaa RE 


AU, 


势能 1 


aE Amm? 
=] 


地 貌 演化 年 份 


田 剑 等 3 的 研究 证 明 沟 密 密 度 和 土壤 侵蚀 强 
度 存在 显著 的 正 相 关 , 随 着 沟 蜜 密度 不 断 增 大 , 土 
壤 侵 蚀 逐 渐 加 强 。 因 此 ,本 文 拟 通过 研究 黄土 小 流 
域 地 貌 系统 的 势能 信息 炳 与 沟 蜜 密 度 之 间 的 相关 
关系 ,间接 验证 势能 信息 炉 与 黄土 侵蚀 之 间 的 关 
系 。 利 用 前 述 方法 计算 出 辛 店 沟 小 流域 5 期 DEM 
的 势能 信息 信 , 然 后 采用 南 学 良 等 的 沟谷 网 络 密度 
二 阶 导 函 数 法 确定 该 小 流域 的 汇流 累积 量 阐 值 为 
12”, 据 此 提取 5 期 DEM 的 沟谷 网 络 并 计算 出 其 沟 
RE ,结果 如 图 5 所 示 。 

从 图 5a 中 可 以 看 出 ,5 期 DEM RRE ARI 
为 负 炉 , 且 总 体 上 旦 不 断 减 小 态势 。 在 整个 地 貌 演 
化 过 程 中 ,势能 信息 炉 最 小 值 为 -1.187, 最 大 值 为 
-1.096。 从 2000 4F F] 2019 E 44 AE IE A a > T 
8.30% 。 与 之 相对 应 , VA) ARF RE RAS E Se 7) N 
的 态势 ,最 小 值 为 1.771 km .km ,最 大 值 为 1.896 
km:km*。 从 2000 年 到 2019 年 沟 蜜 密 度 仅 小 幅 增 
加 7.06%。 这 表明 , 辛 店 沟 小 流域 地 貌 系统 的 演化 
过 程 是 一 个 势能 信息 信 不 断 减 小 的 过 程 ,同时 也 是 
该 系统 内 部 的 沟 密 密度 小 幅 增 大 、 黄 土 侵蚀 小 幅 增 
强 的 过 程 。 对 势能 信息 粹 和 沟 蜜 密 度 进行 多 项 式 
拟 合 ,结果 如 图 5b 所 示 。 从 该 图 可 以 看 出 , 沟 蜜 密 
度 与 势能 信息 炉 之 间 有 具有 良好 的 三 次 多 项 式 函 数 
关系 。 利 用 该 拟 合 公式 ,可 通过 计算 小 流域 地 貌 系 
统 的 势能 信息 箭 来 大 致 估算 其 沟 窄 密度。 这 表明 


(b) BE te SYS SLR EL HR 


1.92 p — ALE SARE 
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Fig. 5 Variation and fitting curve of potential energy information entropy and gully density in Xindiangou small watershed 
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SE WE YAZIN Yat R A BE E I ED A DL A TA tH 
貌 侵蚀 小 幅 增强 的 过 程 ABS MIB WH — 
致 ,进一步 证 实 了 该 结论 的 正确 性 。 
3.2 AR SMM BE 
姜 琼 ”的 研究 表明 ,在 室内 人 工 降 雨 模拟 实验 
条 件 下 , 均 质 黄土 小 流域 的 势能 信息 信 能 够 很 好 地 
指示 黄土 小 流域 的 地 貌 发 育 阶 段 。 为 进一步 验证 
该 结论 是 否 适用 于 野外 非 均 质 多 岩 土 层 黄土 小 流 
域 地 貌 系 统 , 本 文 利用 势能 信息 炉 对 辛 店 沟 黄土 小 
流域 地 貌 系统 的 地 貌 发 育 过 程 进 行人 研究。 首先 , 采 
用 起 伏 比 法 计算 5 期 DEM 的 Strahler 面积 高 程 积 分 
(Hypsometric index, HI) 。 起 伏 比 法 是 Pike 等 "经 
过 数学 公式 推导 得 出 的 估算 HI 的 方法 ,计算 公式 
如 下 : 
HI= (平均 高 程 - 最 小 高 程 )/( 最 大 高 程 z 最 小 高 程 )) 
(6) 
本 文 利用 ArcGIS 软件 提 取 不 同时 期 DEM 的 最 
大 值 .最 小 值 及 平均 值 ,再 利用 式 (6) 求 取 H”™, 计 
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算 结 果 如 表 3 所 示 。 

从 表 3 中 可 以 看 出 ,HI 随 着 地 貌 演化 年 份 的 增 
加 呈现 不 断 减 小 趋势 ,2000 年 HI 最 大 ,为 0.534， 
2019 年 HI 最 小 ,为 0.485,2000 一 2019 年 地 貌 演化 过 
程 中 ,HI 减 小 了 9.18%。5 期 DEM 的 HI 均 位 于 0.4~ 
0.6 之 间 , 因 此 ,从 2000 一 2019 年 辛 店 沟 小 流域 地 貌 
发 育 阶 段 均 为 壮年 期 。 

为 了 更 直观 对 比 室内 小 流域 与 野外 辛 店 沟 小 
Wit AE A hl EP AY HI AA AY SH ASB A 
FEIE , ASCE il T 2h) ita AY ARE a A iS HIS 
化 曲线 ,如 图 6 所 示 。 

姜 琼 中 的 研究 表明 ,在 地 貌 发 育 的 幼年 期 势能 
言 息 炉 不 断 减 小 ,在 壮年 期 势能 信息 炉 不 断 增 大 
(图 6a)。 然 而 , 辛 店 沟 小 流域 的 地 貌 发 育 过 程 表 
HH ,在 其 地 貌 发 育 的 壮年 期 势能 信息 炉 却 不 断 减 
小 ,与 室内 小 流域 势能 信息 信 变 化 不 一 人 至。 为 此 ， 
本 文 结合 戴 维 斯 的 地 貌 侵 蚀 循 环 理论 ”势能 信息 
烂 的 概念 和 室内 小 流域 . 辛 店 沟 小 流域 的 势能 信息 


#3 5 期 DEM $} Strahler 面积 高 程 积分 
Tab.3 Strahler hypsometric index of 5 periods’ DEM 


DEM 数据 地 貌 演化 年 份 平均 高 程 /m 最 小 高 程 /m 最 大 高 程 /m HI 
第 1 期 2000 930 820 1026 0.534 
第 2 期 2007 903 789 1009 0.518 
第 3 期 2009 924 814 1039 0.489 
第 4 期 2015 904 791 1023 0.487 
第 5 期 2019 935 820 1057 0.485 

注 :HI 为 面积 高 程 积分 。 
-7.96 (a) 室内 小 流域 . 0.64 -1.08 (b) 辛 店 沟 小 流域 _ 0.540 
-7.98 PHBE | e -1.09 一 - HA | 0.534 
-8.00 — HI ' -1.10 == HI Ds 
_8.02 0.62 -1.11 l 
图 8.04 0.61 @ -1.12 
m -8.06 0.60 m -1.13 B 
= -8.08 > E Æa T 
ze -8.10 0.59 ge -1.15 
-8.12 0.58 -1.16 
-8.14 0.57 -1.17 
-8.16 : -1.18 
-8.18 0.56 -1.19 
8.20 ! 1 ! 0.55 -1.20 
-10 0 10 20 30 40 50 60 70 80 qu AN ve AS & 6 ab Dd 
Ce 
地 和 貌 演化 天 数 /d 地 和 貌 演化 年 份 


注 :HI 为 面积 高 程 积分 。 下 同 。 


图 6 势能 信息 炉 和 面积 高 程 积分 随地 貌 演化 时 间 的 变化 


Fig. 6 Variations of potential energy information entropy and hypsometric index with geomorphic evolution time 
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灶 变 化 特征 ,对 产生 该 差异 的 原因 进行 分 析 。 假 设 。 为 0.59, 非 常 接近 幼年 期 和 壮年 期 的 分 界 值 0.6。 此 
存在 一 个 地 表 完 全 水 平 的 流域 地 貌 系统 ,在 地 貌 开 ”外 , 辛 店 沟 小 流域 的 势能 信息 炉 在 EF 段 不 断 减 小 ， 
始 发 育 演化 时 ,此 时 处 于 流域 地 貌 发 育 阶段 的 幼年 。 ”根据 前 述 分 析 , 在 该 流域 地 貌 演化 结束 时 , 辛 店 沟 
期 最 初 阶段 ,该 系统 内 部 任意 2 个 点 的 高 程 都 相等 ， ”小 流域 的 势能 信息 炳 应 趋 近 于 0。 因 此 ,本 文大 胆 
所 有 剂 分 单元 格 的 势能 均 属于 同一 个 级 别 。 根 据 “， 推测 辛 店 沟 小 流域 的 势能 信息 焙 可 能 会 在 壮年 期 
APEE MTN RE ,此 时 该 流域 的 势能 信息 烂 应 为 ”和 老年 期 分 界 点 Q 处 发 生 转折 ,在 QL 段 不 断 增加 ， 
0, 即 任意 一 个 单元 格 的 势能 等 级 都 有 100% 的 概率 ”并 最 终 趋 近 于 0。 同 理 ,进一步 推测 室内 小 流域 的 
会 落 入 该 唯一 努 能 级 别 中 。 随 着 地 貌 不 断 演化 , 流 YA EARL AT RES PHA A EH OE Tb 
域内 部 地 表 由 于 受到 降雨 等 侵蚀, 开始 逐步 降低 并 ”发 生 转 折 , 在 TL 有 段 不 断 增加 ,并 最 终 趋 近 于 0。 此 
形成 沟谷 网 络 ,逐步 形成 多 个 不 同 的 势能 等 级 ,使 《外 ,由 于 沿 亚 方 向 辛 店 沟 小 流域 的 势能 信息 箭 不 断 
得 每 一 个 单元 格 的 势能 等 级 的 概率 均 小 于 1, 这 表 增加 ,而 室内 小 流域 沿 RH 方 向 势能 信息 箭 也 是 增 
明 随 着 地 貌 的 不 断 侵蚀 ,势能 信息 和 将 不 断 减 小 。 ”加 的 ,由 此 推测 与 之 对 应 沿 NP 方 向 的 势能 信息 焙 
但 根据 戴 维 斯 的 地 貌 循环 理论 ,该 流域 地 貌 最 终 的 。 二 应 增加 ,从 而 形成 损 点 P, 同 理 , 在 室内 小 流域 曲 
地 名 演化 目标 是 形成 夷 平面 或 准 平原 ,也 就 是 最 终 。 线 上 会 形成 拐点 S。 最 终 ,推测 室内 和 野外 小 流域 
该 流域 地 貌 系统 的 地 表 将 再 次 变 为 水 平 或 近似 水 。 地 有 狐 系 统 在 其 完整 的 地 貌 演化 过 程 中 均 将 形成 类 


平 ,此 时 势能 信息 恼 将 再 次 变 为 0。 本 文 把 上 述 流 似 W 型 的 势能 信息 炳 变化 曲线 ,这 表明 热能 信息 信 


域 地 貌 元 整 演 化 过 程 中 形成 的 势能 信息 入 变化 曲 。 能 够 很 好 地 指示 流域 地 貌 幼年 期 .壮年 期 和 老年 期 


线 称 为 该 流域 的 标准 热能 PARZE H o 根据 的 地 貌 发 育 分 界 点 : 较 好 地 指示 流域 地 貌 发 育 的 
上 述 原理 和 室内 小 流域 及 辛 店 沟 小 流域 势能 信息 阶段 。 
炉 的 变化 情况 ,分 别 推测 出 室内 小 流域 和 野外 小 流 
域 地 貌 演化 的 标准 势能 信息 粹 和 HI 变化 曲线 ,如 图 4 讨论 
7 所 示 。 

从 图 7 中 可 以 看 出 ,室内 小 流域 的 势能 信息 信 黄土 高 原 流域 地 貌 系统 的 地 貌 形态 foe A AZ 
在 CR 段 不 断 减 小 ,在 RH 上 段 不 断 增加 ,有 旦 R 点 的 HI ” 育 特征 由 于 受到 流域 内 部 地 层 岩 性 植被 .降雨 和 


be CoH = ESM EURICEPSEY IC) --- 标准 势能 信息 倘 ( 室 内 小 流域 ) 
K 
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地 貌 演化 阶段 


注 :HI 曲线 上 的 AB 段 和 CD 有 段 分 别 为 室内 小 流域 . 辛 店 沟 小 流域 的 真实 HI; 标 准 势能 信息 箭 曲线 上 的 GRH 段 实 线 .EF 段 实 线 分 别 为 室内 小 流 
域 辛 店 沟 小 流域 的 真实 势能 信息 炉 ;1~Le、Ji~Jis、01~0w 分 别 为 地 貌 演化 阶段 的 幼年 期 .壮年 期 和 老年 期 KK LMN PQS TANHA 
小 流域 及 野外 小 流域 的 标准 势能 信息 焙 变 化 曲线 上 的 起 始 以 及 趋势 变化 转折 点 。 

图 7 室内 和 室外 小 流域 的 标准 势能 信息 粹 和 面积 高 程 积分 变化 


Fig. 7 Variations of the standard potential energy information entropy and hypsometric index of indoor and outdoor watersheds 
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人 类 活动 等 诸多 因素 影响 ,使 其 演化 特征 十 分 复 
杂 。 为 此 ,本 文 充 分 考虑 辛 店 沟 小 流域 的 地 层 岩 性 
和 岩 土 层 密度 等 因素 ,通过 构建 辛 店 沟 小 流域 地 貌 
系统 及 其 势能 信息 炉 的 数学 模型 ,从 流域 地 貌 系统 
的 势能 信息 炉 这 一 独特 视角 对 黄土 小 流域 的 地 貌 
演化 的 综合 特征 进行 探索 。 上 述 研究 结果 表明 ,本 
文 构建 的 辛 店 沟 多 兰 土 层 流域 地 貌 系统 及 其 势能 
Fi WBE A Ar OT SE JE YA 7) it Sak AE IT RUBER W o 
野外 自然 条 件 下 辛 店 沟 小 流域 的 地 貌 演化 过 程 是 
AREA BSR EE, Aaa eS > TES AY 
人 工 降 雨 模拟 试验 条 件 下 的 均 质 黄土 地 貌 演化 过 
程 研 究 中 得 出 的 结论 一 臻 。 这 说 明 势能 信息 信 的 
炉 减 过 程 能 够 有 效 反 映 黄土 小 流域 地 貌 不 断 侵蚀 
的 过 程 。 此 外 ,本 文 提出 的 利用 标准 势能 信息 炉 曲 
线 对 流域 地 貌 发 育 阶段 进行 划分 的 方法 ,尽管 还 有 
待 于 后 续 开 展 更 深入 研究 进行 验证 ,为 流域 地 貌 发 
育 阶 段 的 精细 划分 提供 了 一 条 新 的 途径 ,非常 适合 
处 于 不 同 发 育 阶段 的 流域 地 貌 对 比 研 究 。 未 来 我 
们 将 选取 更 多 的 处 于 幼年 期 .壮年 期 和 老年 期 等 不 
同 发 育 阶 段 的 黄土 小 流域 进行 深入 人 研究 ,对 标准 势 
能 炉 曲 线形 态 、 曲 线 拐点 的 地 貌 演化 意义 等 进行 
验证 。 


5 结论 


(1) 本 文 基于 前 人 研究 及 辛 店 沟 小 流域 野外 实 
地 考察 获取 的 数据 ,成 功 构 建 了 野外 多 岩 士 层 黄土 
小 流域 地 貌 系统 的 概念 模型 及 其 势能 信息 炉 的 数 
学 模型 。 该 模型 能 够 有 效 对 辛 店 沟 小 流域 进行 数 
(EBS, AT A IY i A RA fa E A 
学 模型 的 扩展 ,可 在 非 均 质 多 岩 土 层 流 域 地 貌 系统 
中 使 用 ,具有 更 广泛 适用 性 。 

(2) 研究 结果 表明 ,以 黄土 侵蚀 作用 为 主 的 辛 
店 沟 小 流域 从 2000 一 2019 年 的 地 貌 演化 过 程 是 其 
SAN BE fs E I AO RA ak ash EE AAT ETS Sh AN WB e h P 

(3) SEI YA 7] iit Sak AY AAEN AA FI FE AW 8 AN 
该 流域 地 貌 发 育 的 壮年 期 和 老年 期 分 界 点 ,室内 小 
Trt Sa ES ARENE E EA BC I-A A E 
分 界 点 。 为 此 ,本 文 推 测 : 黄 土 流域 地 貌 系统 在 其 
完整 的 地 貌 演化 过 程 中 将 形成 类 似 W 型 的 标准 势 
能 信息 箭 变化 曲线 ,该 曲线 能 较 好 指示 流域 地 貌 幼 
年 期 壮年 期 和 老年 期 的 地 貌 发 育 分 界 点 ,可 为 流 
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域 地 貌 发 育 阶段 的 精细 划分 提供 一 条 新 途径 。 
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Geomorphic evolution characteristics of small loess watersheds based on 


potential energy information entropy 


ZHAO Weidong, WANG Shugin, TIAN Jian, JI Bin, MA Lei 


(School of Resources and Environmental Engineering, Hefei University of Technology, Hefei 230009, Anhui, China) 


Abstract: Research on the evolution characteristics of loess erosion and the development of watershed 
geomorphic systems in the Loess Plateau of China is a hotspot in the study of loess landforms, and many related 
scientific issues still need to be studied. Most previous studies have focused on the erosion and development 
characteristics of watershed geomorphology, but there is a lack of research on indepth analysis of the geomorphic 
evolution characteristics of loess small watersheds with multiple rock and soil layers in the field from the 
perspective of watershed geomorphology system and its potential energy information entropy (PEIE). Therefore, 
based on the viewpoint and method of system theory, mathematical models of the geomorphic system and its 
PEIE of a small loess watershed with multiple rock and soil layers were built. Taking the Xindiangou small 
watershed at Suide County, Shaanxi Province, China as the research sample area, the digital elevation model data 
of 5 periods from 2000 to 2019 were used to investigate the entropy change law of the PEIE and the 
characteristics of landform evolution in the watershed. The results show that from 2000 to 2019, the geomorphic 
evolution process dominated by loess erosion in the Xindiangou small watershed is the process of entropy 
reduction of its PEIE and continuous erosion of the loess landform. Combined with the entropy change of PEIE 
in the Xindiangou small watershed and the indoor small watershed, it is speculated that the watershed 
geomorphic system will form a W-shaped PEIE change curve in its complete geomorphic evolution process, and 
its PEIE can better indicate the watershed geomorphic development demarcation points in its infancy, maturity, 
and old age. The research results also confirm the effectiveness of the conceptual model of the loess watershed 
geomorphic system and its mathematical model of PEIE in the numerical simulation of the Xindiangou small 
watershed. This model is an extension of the existing mathematical model of PEIE of watershed geomorphic 
systems with homogeneous single loess layers, which can be used in heterogeneous multiple rock and soil layer 
loess watershed geomorphic systems and has wider applicability. The research results provide ideas for further 
research on the geomorphic formation mechanism and evolution law of the Xindiangou small watershed in the 
future and guide soil and water conservation, ecological restoration, and regional sustainable development of the 
Xindiangou small watershed, which has important theoretical significance and good application prospects. 

Key words: digital elevation model (DEM); digital terrain analysis; potential energy information entropy; wa- 


tershed geomorphic system; Xindiangou small watershed 


